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There is a number of explicit kinetic energy density functionals for non-interacting electron sys- 
tems that are obtained in terms of the electron density and its derivatives. These semilocal func- 
tionals have been widely used in the literature. In this work we present a comparative study of the 
kinetic energy density of these semilocal functionals, stressing the importance of the local behavior 
to assess the quality of the functionals. We propose a quality factor that measures the local dif- 
ferences between the usual orbital-based kinetic energy density distributions and the approximated 
ones, allowing to ensure if the good results obtained for the total kinetic energies with these semilo- 
cal functionals are due to their correct local performance or to error cancellations. We have also 
included contributions coming from the laplacian of the electron density to work with an infinite 
set of kinetic energy densities. For all the functionals but one we have found that their success in 
the evaluation of the total kinetic energy are due to global error cancellations, whereas the local be- 
havior of their kinetic energy density becomes worse than that corresponding to the Thomas-Fermi 
functional. 
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I. INTRODUCTION 

Density Functional Theory (DFT) has nowadays a 
privileged position within the current methods for cal- 
culating the electron structure. The theorems of Hohen- 
berg and Kohn^ show that the ground state of a system 
of nuclei and electrons can be fully described in terms 
of its electron density. In fact, the total energy E[n] of 
the electrons can be written as a functional of the den- 
sity and the variational minimization of E[n\ yields the 
ground state electron density and its total energy. 

The usual procedure within the DFT solves the Kohn- 
Sham^ (KS) equations for N orbitals, where N is the 
number of electrons of the system. This KS method al- 
lows the partition of the total energy functional in differ- 
ent pieces with distinct physical meanings 

E[n] = Ts[n] + V[n] + J[n] + E,,[nl (1) 

where E[n\ is the total energy, V[n\ is the energy of the 
electron density allocated in the electric field generated 
by the nuclei, J[n] is the classical repulsion of the electron 
density (also called Hartree energy), Ts[n] is the kinetic 
energy of a noninteracting system that yields the same 
electron density of the interacting one and i<^a;c[n] is a 
density functional that takes into account the exchange 
and the correlation (XC) energies. The major part of 
the kinetic energy, T5[n], is then expressed exactly by 
means of the KS one-electron orbitals and the small part 
(T[n] — ^^[n]) is included in the XC functional, Exc[n], 
which is the only part of the energy to be approximated. 
The exact ground-state electron density n (r) is the sum 
of the densities of the KS orbitals, the kinetic energy 
functional is evaluated exactly through the KS orbitals,^ 
and the minimization of the total energy functional ([T]) 
becomes the resolution of a set of coupled Schrodinger- 



like equations for the KS orbitals. 

The computational cost of solving the KS equations 
can be avoided by using Ts as a functional depending 
explicitly on the electron density, instead of construct- 
ing the KS orbitals and evaluating the kinetic energy in 
terms of them. In that case, the computational cost of 
this orbital-free procedure would scale with the number 
of electrons (a linear scaling method) and can offer much 
faster computational times than the KS calculations, al- 
lowing to deal with systems that involve several hundred 
thousands of atoms or more. The ground state of these 
systems is then calculated through an Euler-Lagrange 
variational minimization. Moreover, the use of an ap- 
proximate orbital-free kinetic energy functional, instead 
of the KS method, makes easier the evaluation of the 
forces and reduces the computational cost of complex 
calculations in first-principle molecular dynamics'^. But 
one of the issues that must be solved about orbital-free 
kinetic energy functionals is the stability of their numeri- 
cal solutions, because some functionals have been proved 
to be linearly stable but nonlinearly unstable (see, e. g., 
Ref. 4) and the solution obtained can be meaningless. 

As a consequence, the construction of functionals de- 
pending explicitly on the density n (r) (without any ref- 
erence to the wave functions) has an undoubted for- 
mal interest and an important practical side. Being 
the analytic dependence of Ts on the electron density 
a rather academic problem for a long time (see re- 
views m Refs. H and H), the increasingly rapid devel- 
opment of computational chemistry has turn it an in- 
teresting topic, from the proposals of new kinetic en- 
ergy functionalsii^^iaOiiisasii^ie^isas to the study 
of the kinetic energy density^"'^^ and the application to 
simple systems i^^i^i^i 

In order to construct accurate explicit Ts [n] function- 
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als we need to get not only good total energies but also 
the correct density profiles for the ground state of the 
total energy functional (Eq. [1]) after a fully variational 
minimization of the energy. For that reason, although 
common tests to determine the quality of a kinetic en- 
ergy functional only calculate the energies and density 
profiles using good densities (i.e., those obtained with ac- 
curate methods such as the Hartree-Fock or KS ones) we 
are interested in studying the quality of the kinetic func- 
tional beyond this limitation. In particular, we think the 
properties of the kinetic energy density deserve a study 
by itself, paving the way for understanding when and why 
functionals are able or not to describe the characteristic 
quantum properties, like the presence of structure in the 
density profiles. 

This paper presents an extensive study of the kinetic 
energy densities of a number of semilocal kinetic func- 
tionals thanks to the proposal of a quality factor that 
helps to determine how far from a valid kinetic energy 
density are the approximated ones. 



II. THE KINETIC ENERGY DENSITY 

The kinetic energy density (KED) can be defined as 
any function ts{r) that integrates to the exact total ki- 
netic energy, 

Ts[n] = Jdr tsir). (2) 

It is clear that such definition does not determine the ki- 
netic energy density uniquely: any function - with the 
appropriate scaling properties - that integrates to zero 
in the whole space can be added to any KED to yield 
another KED, and an infinite number of new valid KEDs 
can be defined by multiplying this function by any coef- 
ficient. The non uniqueness of the kinetic energy density 
has been studied in the literaturci^i^ 

Within the KS method, an orbital-based KED can be 
obtained in terms of the KS orbitals. Two main defini- 
tions are commonly used in the literature. The first def- 
inition has the advantage of being positive everywhere, 
because it is calculated through the squared gradient of 
the orbitals of an iV-electron system (atomic unit will be 
used in this paper): 

1 ^ 

4{r) = -J2\\7Mr)\\ (3) 

1=1 

where (piir) is the ith KS orbital. A second definition can 
be obtained using the kinetic energy operator in the way 
it appears in the KS equations, yielding a non positive 
definite function: 

1 ^ 

4'(r) = -2E'^nr)V^'/'.(r). (4) 



Functions tg and tg^ are related through the laplacian of 
the electron density, 

to(r) - V2n(r), (5) 

as 

4(r)-4^(r) = ito(r). (6) 

They have different local properties, but both of them are 
valid definitions for KED because they integrate to the 
same total kinetic energy since the integral of to(r) over 
the whole space is zero. In finite systems the electron 
density decays exponentially and the divergence theorem 
can be used to transform the integral over the space into a 
surface integral, (irV^7i(r) = /^(iSVn(r). The gradi- 
ent of the electron density decays faster than the growth 
of the surface (which is proportional to r^) and the in- 
tegral is zero. For extended systems, but periodic, the 
last integral is extended over the surface of the unit cell 
of the periodic system; the gradients in opposites sides 
of the cell cancel one each other, and the surface integral 
has also a zero value. 

We must remark that, besides the interest of the 
study of the intrinsic properties of the KED, some new 
XC density functionals have been proposed using the 
KED^>26^ 



III. SIMPLE KINETIC ENERGY DENSITY 
FUNCTIONALS 

A. The First Functional Approximations. 

The first explicit kinetic energy density func- 
tional formulated was the Thomas-Fermi (TF) 
approximation, where each point of an inhomoge- 
neous electron system with density n(r) contributes to 
the kinetic energy as any point in an homogeneous one 
with the same electron density, no — n,{r). We then get 
the TF functional 

TTF[n{r)] = J drCTPn'/^r), (7) 

where Ctf — ^(Stt^)^/^ is the TF constant. By con- 
struction, this functional is exact for homogeneous sys- 
tems and a good approximation for systems close to the 
free electron gasi^i^ But this functional gives, in gen- 
eral, poor results when applied to atoms or molecules, 
and the density profiles obtained in a variational mini- 
mization show no quantum effects. For atoms, in fact, 
no shell structure is obtained and the errors in the to- 
tal kinetic energy are about 10%. This relative error is 
also obtained when the functional is applied using good 
densities. 

Another explicit kinetic energy density functional is 
the von Weizsacker (vW) Tuvi/in] functional, constructed 
to be exact for one or two electrons in the same spacial 
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state. ^ The most usual form in which functional is found 
in the literature is 



Tvw[n] = 



|Vn(r) 
n(r) 



(8) 



but can also rewritten in other ways that yield the same 
total kinetic energy. Being the functional exact for sin- 
gle orbital systems, it gives the correct KED when the 
contribution of a given orbital to the electron density is 
much larger than all the other orbitals. For the same rea- 
son, for those regions in space that are very close or very 
far away from the positions of the nuclei, this functional 
also gives the correct KED. Finally, Ti,iy[n] is exact in 
regions with large variations of the electron density. But 
when applied to general systems this functional yields 
large errors. So, the variational minimization for atoms 
gives no shell structures and the relative errors increase a 
lot when the number of electrons grows, getting energies 
of a different order of magnitude than the exact ones for 
atoms with a large number of electrons. 

A very interesting functional can be obtained as a com- 
bination of the TF and vW functionals, the so called 
second order gradient expansion approximation (GEA2), 
constructed as>^ 



GEA2! 



Ctf J n(r)5/3rfr+l J 
Ctf j dvn{rf/ 



|Vn(r) 
n(r) 



-dr (9) 



TICtf 



Vn(r) 



24/3 (r 



This functional gives the correct kinetic energies for sys- 
tems with slow varying electron densities. Density pro- 
files obtained in the variational minimization of the to- 
tal energy with the GEA2 approximation show the same 
pathologies as the previous functionals, but when applied 
using good electron densities the error in the total energy 
is only about 1% for a number of very different systems. 
This error is clearly too big for chemical precision, but is 
quite small for common orbital-free kinetic energy den- 
sity functionals 



B. The Generalized Gradient Approximations. 

In this paper we make a comprehensive study of the 
kinetic energy density functionals that can be expressed 
within the Generalized Gradient Approximation (GGA), 
which allows a general form for those semilocal function- 
als that only depend on the electron density and its gradi- 
ent. Any semilocal functional of this kind can be written 
as 



GGA\ 



dYtTF{r)F,,,h{s{Y)), (10) 



where Itf (r) = CTF'T-^^'^(r) is the KED corresponding 
to the TF functional and Fe„/i(s(r)) is the so-called en- 



hancement factor, that depends on the adimensional vari- 
able 



s{v) 



Vn(r)| 
n4/3(r) ■ 



The quantity s(r) is called the reduced density gradient 
and has a clear physical interpretation, because it con- 
trols the speed of the variation of the electron density. 
Large values of s(r) correspond to fast variations in the 
electron density and small values to slow ones; a zero 
value indicates a region of the space where the electron 
density has no variation. 

The mathematical form of the enhancement factor de- 
termines the functional and authors have proposed many 
different forms for the enhancement factor. In this pa- 
per, the following set of semilocal functionals have been 
selected. 



1. Thomas-Fermi functional (TF)3>^^ When the 
enhancement factor has a constant value of 1, the TF 
functional is recovered, 

Ftf{s{y)) = 1. 

2. Second order gradient expansion approxima- 
tion (GEA2).'^'^ This approximation is a particular case 
of the GGA approximations, obtained with the expres- 



FgeA2{s{v)) 



1 



72Ctf 



There are some functionals constructed as a linear 
combination of the TF functional and the vW functional. 
We cite now some of them. 

3. Thomas Fermi + ^ von Weizsacker (TF5W). 
In this case the contribution of the vW functional is 
weighted by a prefactor l/5^^^^^2^^^^^^ 



F 



TF5W 



(sir)) = 



1 



1 



AOCtf 



sir? 



4. Thomas Fermi + von Weizsacker (TFvW)j^ 

The combination involving both the full Thomas-Fermi 
and von Weizsacker functionals is theoretically interest- 
ing but in practical applications largely overestimate the 
kinetic energy. 



Fr 



TFvW 



(sir)) 



1 



1 



8Ctf 



sir? 



5. Thomas Fermi + | von Weizsacker 
(TF9W). A modification of GEA2 can be written with 
the introduction of a parameter b that modifies the con- 
tribution of the gradient correction. Several values of b 
have been used by different authors. For this work we 
have chosen a value of 6 = 1.0674^ 



F, 



TF9W 



(sir)) = 



1 



72Ctf 



sir? 
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6. A^-dependent Thomas Fermi functional (TF- 
N)i^ Several functionals have been developed as modifi- 
cations of tlie TF functional by including a prefactor that 
depends on the number of electrons N. For this work we 
have chosen: 



0.313 0.187 



FTF-Nisir)) = 1 + 



iVi/3 Ar2/3 



that usually provides better values for the kinetic energies 
than those obtained with the TF functional. 

7. Pearson functional (Pear).^^ This functional, a 
modification of the GEA2 approximation that follows the 
idea presented by Pearson and Gordon^ is constructed 
in such a way that the gradient correction takes only into 
account the regions of the space where the density varies 
slowly. This can be done with the introduction of a cut- 
off that depends on the value of the reduced gradient 
s. A sharp cutoff usually introduces numerical problems 
and a functional with a smooth cutofF'^'^ is proposed with 
the enhancement factor: 



Fpear [n] = j '^^ 



1 



1 



1 + 72Ctf 



where C is a parameter that Pearson and Gordon fixed as 
C = 1, the optimum value for fitting to atomic Hartree- 
Fock kinetic energies. This functional form prevents the 
contributions of the gradient expansions coming from re- 
gions where fast variations of the electron density occur. 

8. DePristo Kress functional (DK)4^ These au- 
thors introduce a mathematical form based on a Pade- 
type approximation that is able to recover different de- 
sirable limits: 



FDK{s{r)) = 



a^x^ + a2x'^ + aix + 1 



bsX'^ + b2x'^ + bix + 1 



where 

72Ctf 

This functional is equivalent to GEA2 for slowly vary- 
ing densities - small values of s - and equivalent to the 
von Weizsacker functional for fast varying densities - 
large values of s The parameters of the functional 
(ai = 0.95, 02 = 14.28111, = -19.5762, hi = -0.05, 
62 = 9.99802, bs = 2.96085) were obtained by fitting the 
results of the functional to atomic Hartree-Fock kinetic 
energies: 

9. Lee Lee Parr functional (LLP)."'^ These au- 
thors have used the same mathematical form of an en- 
hancement factor for the exchange functional for con- 
structing the kinetic energy density functional 



FLLpisir)) 



6[s(r) 



1 + c [s(r)] arcsin(s(r)) 



where the parameters have the values b = 0.0044188 and 
c = 0.0253. 

10. Ou-Yang Levy 1 functional (OLl)4^ Fol- 
lowing the nonuniform coordinate scaling requirement 
for the kinetic energy density functional;^ these authors 
suggest these two new functional forms: 

FoLiisir)) = 1 + [s{r)f + d [.(r)] , (12) 



TF 



where d ^ 0.00187, and 

11. Ou-Yang - Levy 2 functional (OL2)4^ 



(.(r)) = l + 



1 



TICtf 



D[s{v)] 



1-H25/3 [s(r)]' 



with D = 0.0245. 

12. Thakkar functional (Th.ak)M In a review of 
many semilocal functionals, Thakkar proposed a new one 
as a conjoint of the mathematical forms of the most suc- 
cessful functionals: 



FThak{s{v)) = 1 + 



bW)Y 



D [s{v)] 



1 + c [s(r)] arcsin[s(r)] 1 + 2^5/3 [s(r)] ' 



with the parameters: b = 0.0055, c — 0.0253 and D = 
0.072. 

As was commented before, Lee, Lee and Parr argued 
that the enhancement factors used for the GGA XC en- 
ergy functionals can also be used for the development of 
kinetic energy functionals. The next five functionals, for- 
mulated following that approach, were proposed by Lacks 
and Gordon.'^*' 

13. Becke 86A functional (B86A)^ 



FB86A{s{r)) = 1 + 0.0039: 



1 + 0.004 [s{r)Y 
14. Becke 86B functional (B86B)^ 



FB&6Bis{r)) = 1 + 0.00403. 



Hr)Y 



1 + 0.007 [s{r)Y] 
15. DePristo Kress 87 functional (DK87)^ 



DK87 



(sir)) = 1+ 



324(187r4)^/^ 



Kr)] 



2 1 + 0.861504s(r) 



1 + 0.044286 [sir)Y 

16. Perdew-Wang 86 functional (PW86)^ 

Fpwseisir)) = 1+1.296 [s(r)] Vl4 [s(r)]4+0.2 [s{r)f/'' 

17. Perdew Wang 91 functional (PW91)^ 

Fpw9i{s{r)) = 



1 + ais(r) argsinh (6 [s(r)]) 



0-2 — 036 



-100[s(r)]^ 



1 + as(r) argsinh (6 [s(r)]) + 04 [s(r)]^ 
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with ai = 0.19645, 03 = 0.2747, 03 = 0.1508, 04 = 0.004 
and b = 7.7956. 

18. Lacks Gordon functional (LG94). In 

their extensive study of the kinetic energy density 
functionals)^ these authors made their own contribution 



with the formulation of a new one, based on a previous 
exchange functional of the same authors developed by fit- 
ting exchange energies of atoms and ions, as well as the 
correct behavior of the exchange energy for small s(r): 



^LG94(s(r)) = 



1 + 02 [s(r) 



■ a4 [s{r)f + ae [s{r)f + as [sir)f + aio [s(r)]^ 



- ai2 [s 



1 + 10-8 [s{r)Y 



r 



29.790, ae = 
= 55.944 and 



with aa = (iQ-^ + 0.1234) /0.024974, 04 = 
22.417, as = 12.119, aio = 1570.1, 012 
b = 0.024974. 

19. von Weizsacker functional (vW).'^^ The vW 

functional can also be written as a semilocal functional 
if we choose the enhancement factor: 



Fywisir)) 



1 



8G 



TF 



-s{vf. 



20. Acharya et al. functional (ABSP).^^ 

Acharya, Bartolotti, Sears, and Parr proposed a func- 
tional with the full vW functional, introducing a TF 
contribution weighted by a prefactor depending on the 
electron number N: 



Fabsp[iA = 



1 



8C' 



TF 



-s(r) 



1.412 



21. Gazquez Robles functional (GR).^^ Follow- 
ing the same spirit of the previous functional, Gazquez 
and Robles developed another one, with a more compli- 
cated form for the weight of the TF functional: 



FgrH 



1 



8Ctf 



1 



2 

N 



1 



1.303 

iVl/3 



0.029 \ 

iV2/3 ) ■ 

(13) 



IV. MEASUREMENT OF THE DISTANCE 
BETWEEN TWO KINETIC ENERGY 
DENSITIES. 

As the KED can be evaluated by using the KS orbitals, 
we may think that every approximate functional yields 
an approximation to the KED by simply defining it as 
the integrand t^""'^(r) that appears in the expression of 
the functional: 



dr t 



June 



(r). 



(14) 



As the main aim of this paper is to know, in a quan- 
titative way, how the semilocal functional are able to 
reproduce the KED, we need to define a measure of the 



closeness between the orbital-based and approximated 
orbital-free KED distributions, i.e. we must find an ef- 
ficient and reasonable measure of the distances between 
distributions. Some distances between probability distri- 
butions have been defined, most of them to be used for 
normalized and positive definite distributions. The elec- 
tron density is normalized and positive, but the KEDs 
are no normalized and they are not, in general, positive 
definite (note that the second definition of the KS KED 
has always negative values in some regions of the space). 
For these reasons those distances are not appropriate in 
our case. But the absolute value norm, being the ab- 
solute value of the difference between two distributions 
integrated over the whole space, is a measure of the ac- 
cumulated local error of the KED, 



dr 



tsir) 



.June 



(r) 



(15) 



where ts'(r) is a valid distribution of KED and tg""^(r) 
an approximated one. This definition satisfies the three 
desirable requirements for a distance: it is positive defi- 
nite, is symmetric and verifies the triangular inequality. 

The defined measure e depends on the size of the sys- 
tem, because the total kinetic energy has a nonlinear de- 
pendence with the number of electrons (the kinetic en- 
ergy of neon is more than two hundred times the kinetic 
energy of hydrogen) . In order to get a measure that does 
not depend so much on the size of the system we divide e 
by the exact total kinetic energy. In this way, we propose 
a quality factor a for a KED as 



(16) 



Ts\n\ 



With this definition, the bigger the difference between the 
distributions is the larger the value of a becomes, and tr 
is zero for two identical KEDs. Note that the value of 
a can be interpreted as the amount of the approximated 
KED that is misplaced when compared to the distribu- 
tion of the orbital-based KED. With this factor tr we 
are going to test whether the good results obtained for 
the total kinetic energies with the GGA functionals are 
due to their correct local KED or to error cancellations 
among different regions of the space. 
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V. ELECTRON DENSITIES FOR LIGHT 
ATOMS 

In order to clarify the quality of the KED each ap- 
proximate functional yields, we want to use good enough 
atomic electron densities to evaluate our semilocal func- 
tional. 

We initially tried commercial software that performs 
KS calculations and outputs the orbitals used; in that 
way the related electron density, the KED and the total 
kinetic energy can be evaluated. But we decided to avoid 
the use of these codes for several reasons. Firstly, the be- 
havior of any XC functional leads to a bad decay of the 
density in the outer region of the atom, reflecting some 
problems in the quality of the KS orbitals, their density 
and in the KED. More importantly, when using a basis 
set of gaussians, several gaussians are generally needed 
to approximate one Slater-like KS orbital. As a result, 
changes in the slope of the KS orbitals appear; that could 
not noticeable by simple inspection but, when the care- 
ful numerical evaluation of the laplacian of the density is 
done, spurious structure arises. So, when the laplacian 
of the density is numerically evaluated one must be cau- 
tious with the use of a gaussian basis set. Furthermore, 
gaussian basis set are unable to give a good description 
of both the cusps of the density in the positions of the 
nuclei and the decay of the orbitals far away from them. 

For that reasons, we have decided to use Slater orbitals 
for the light atoms. They only approximately describe 
the KS orbitals but all the aforementioned pathologies 
will not appear (there are no spurious oscillations, no 
cut-offs and correct cusp conditions and decay of the or- 
bitals can be achieved). There are several sets of values 
of exponents for the Slater orbitals, with only small dif- 
ferences between them; Table |T] presents our choice of 
parameters^Si^ for the ten first atoms. 



VI. TOTAL KINETIC ENERGIES 

We have performed total energy calculations for twenty 
one functionals, using the densities corresponding to the 
Slater orbitals with parameters given in Table [D 

The relative errors for the total kinetic energies are 
presented in Table [III We also show the average of the 
absolute values of the errors for every functional. To 
add more information, the dispersion of data is displayed 
under the label "Range" , where we give the statistical 
range, i.e., the length of the interval which includes all 
the former data. When the value of the range is larger 
than the average we can think the average could be mean- 
ingless. All the tables of this paper will be presented this 
way (in each case, for the relative error, for the values of 
cr, etc.). 

As we can see in Table we get errors bigger than 

8% for the Ttf, TtF1/5i;IV: TtFvW, Tpear, T^w, Tabsp 

functionals, about 4.5% for the Tgb. one and smaller than 
2% for all the other functionals. 



VII. THE KINETIC ENERGY DENSITY AND 
THE ROLE OF THE LAPLACIAN OF THE 
ELECTRON DENSITY 

The nonuniqueness of the orbital-based definition of 
the KED have been pointed out by many authors. The 
first definition ^^(r) is widely used in the theory Atoms 
in Molecules^^ of Bader, that recently claimed the useful- 
ness of the topological analysis of that definition^. The 
asymptotic behavior of the of the KED seems to support 
the use of the second definition t^g (r) But "classical" 
properties are recovered more appropriately with a mix- 
ture of both definitions, in particular the mean average 
of them (see, e. g., Refs. [63 and [63|) . We feel the use of 
one or another definition is still an open question. 

When we compare the approximated orbital-free KED 
with t^gir) very large values for the quality factor <j are 
obtained, about 0.6 for the TF functional and for almost 
all the other GGA functionals. But, as commented, the 
comparison of the GGA functionals with any single defi- 
nition of KED is only a choice. It is possible to make com- 
parisons with an infinite set of KEDs. Indeed, the differ- 
ence between t^s{Y) and t^g {y) is one fourth the laplacian 
of the electron density, an archetypical function related 
with the electron system that integrates to zero over the 
whole space and having the adequate scaling properties 
to be a kinetic energy density. So, an infinite set of valid 
KEDs can be obtained through (see, e.g., Ref. 20) 

t^(r)=4(r)-Haio(r), (17) 

where the ^^(r) is constructed as ^^(i") plus the laplacian 
of the electron density io(r) — V^n(r), multiplied by a 
prefactor a that can have any real value. The value of 
a — Q yields ^^(r), whereas the value of a — — i gives 
tjJ{Y). On the other hand, the arithmetic mean of the 
first and second definition is recovered with a = 
and we have all the intermediate KEDs by continuously 
varying a. 

We want to test every functional in the more adequate 
conditions for itself. For that reason, we have compared 
its approximate KED with the ^^(r) that is its closest, 
choosing among all possible values of the parameter a. 
To do that we minimize the value of a when varying 
the parameter a. To obtain the value of a that makes 
the best fit of the distribution t^{v) to the approximated 
KED we have minimized the value of a using a golden 
search algorithm^^ After the minimization process, each 
functional has a given value of a that yields the lower 
value of cr, being the closest KED that t^{v) constructed 
by using this value of the parameter a in Eq. (|17p . 

This methodology is used for all the functionals. The 
best values of a are shown in Tabic Hill We see that the 
values of a are almost constant and do not depend so 
much on the number of electrons for almost all of the 
functionals. For all functionals but those with a full vW 
term, the best fits to (r) are almost equidistant to both 
tg(r) and tg^ {r), although a little bit closer to the second 
one. The vW functional is a special case: the values of a 
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H 


He 


Li Be 


B 


C 


N 





F 


Ne 


Is 


1 


1.6875 


2.6906 3.6848 


4.6795 


5.6727 


6.665 


7.6579 


8.6501 


9.6421 


2s 






1.2792 1.912 


2.5762 


3.2166 


3.8474 


4.4916 


5.1276 


5.7584 


2p 








2.4214 


3.1358 


3.834 


4.4532 


5.1000 


5.7584 



TABLE L Values of the effective charge, Z^ff. 





H 


He 


Li 


Be 


B 


C 


N 


O 


F 


Ne 


Ave 


Range 


TF 


-0.082 


-0.082 


-0.085 


-0.091 


-0.097 


-0.095 


-0.085 


-0.090 


-0.087 


-0.077 





087 





020 


GEA2 


0.029 


0.029 


0.023 


0.013 


0.003 


0.000 


0.006 


-0.003 


-0.003 


0.004 





Oil 





032 


TF5W 


0.118 


0.118 


0.109 


0.097 


0.083 


0.077 


0.079 


0.067 


0.064 


0.068 





088 





054 


TFvW 


0.918 


0.918 


0.887 


0.846 


0.800 


0.762 


0.736 


0.695 


0.667 


0.648 





788 





269 


TF9W 


0.036 


0.036 


0.030 


0.020 


0.010 


0.007 


0.013 


0.003 


0.003 


0.009 





017 





034 


TF-N 


0.034 


0.038 


0.031 


0.021 


0.011 


0.010 


0.018 


0.010 


0.011 


0.020 





020 





028 


Pear 


-0.070 


-0.070 


-0.074 


-0.079 


-0.085 


-0.083 


-0.073 


-0.078 


-0.075 


-0.064 





075 





021 


DK 


0.031 


0.031 


0.016 


-0.003 


-0.011 


-0.010 


-0.001 


-0.007 


-0.005 


0.005 





012 





042 


T T Ti 

LLF 


0.023 


0.023 


0.018 


0.011 


0.002 


0.001 


0.008 


0.000 


0.001 


0.008 





010 





022 


OLl 


0.039 


0.039 


0.032 


0.023 


0.012 


0.010 


0.016 


0.006 


0.006 


0.013 





020 





033 


OL2 


0.036 


0.036 


0.029 


0.020 


0.009 


0.007 


0.013 


0.004 


0.004 


0.011 





017 





032 


Thak 


0.026 


0.026 


0.022 


0.015 


0.006 


0.005 


0.011 


0.004 


0.004 


0.011 





013 





022 


B86A 


0.023 


0.023 


0.019 


0.013 


0.004 


0.003 


0.010 


0.002 


0.003 


0.010 





Oil 





021 


B86B 


0.049 


0.049 


0.045 


0.038 


0.028 


0.026 


0.033 


0.025 


0.025 


0.032 





035 





025 


DK87 


0.030 


0.030 


0.026 


0.019 


0.010 


0.009 


0.015 


0.008 


0.008 


0.015 





017 





022 


PW86 


0.020 


0.020 


0.016 


0.009 


0.001 


0.000 


0.007 


-0.001 


0.000 


0.008 





008 





021 


PW91 


0.026 


0.026 


0.021 


0.014 


0.005 


0.003 


0.009 


0.001 


0.001 


0.008 





Oil 





025 


LG94 


0.024 


0.024 


0.020 


0.013 


0.005 


0.004 


0.010 


0.003 


0.003 


0.011 





012 





021 


vW 


0.000 


0.000 


-0.027 


-0.064 


-0.104 


-0.143 


-0.179 


-0.215 


-0.247 


-0.275 





125 





275 


ABSP 


-0.378 


-0.111 


-0.008 


0.037 


0.054 


0.059 


0.061 


0.052 


0.047 


0.043 





085 





439 


GR 


0.252 


0.000 


0.006 


0.023 


0.031 


0.033 


0.035 


0.028 


0.023 


0.022 





045 





252 



TABLE 11: Relative errors in the total kinetic energy for the semilocal functionals. The average is made over the absolute 
relative errors. 



are zero for the four first atoms and very close to zero for 
all the others, reflecting the fact that T^w in the form 
usually found in literature is directly related to tg (r) and 
its approximate KED is always closer to that definition 
than to any other one. 

In Table IIVI we present the corresponding values of 
a. As an average, the TF functional puts a 16.5% of 
the KED misplaced from those regions where the KED 
that better fits TF is located. Thinking that the GGA 
functionals are corrections to the TF functional, one can 
expect that they will improve not only the total energies 
(as they do) but also the local behavior of the kinetic 
functional. GEA2 do the same for 18.8%, and a careful 
exam of the results show an unexpected result: all the 
GGA functionals, despite their improvement in the eval- 
uation of the total energies, yield larger values for the 
quality factor a that the TF ones. It seems that, within 
the GGA scheme, all functionals but those with a full vW 
term improves the TF results giving total kinetic energies 
within 2% of the exact one, but they place the additional 
KED (i.e. the KED not included in the TF functional) in 
wrong regions of the space. We conclude that the GGA 
functionals improve the results for the TF kinetic energy 
by global error cancellations in the evaluation of the total 
kinetic energies, while the local behavior of their KEDs 
becomes worse than the TF one. The only exception 



to the previous results is the Pearson functional, con- 
structed with a different philosophy, which gives slightly 
better values for a and for the total kinetic energies. 

For the purpose of gaining insight on the origin of the 
previous results, we have also divided the contribution 
to the KED coming from each orbital. Even with no 
clear physical justification, we have used a sum of the 
laplacian for each orbital (f>i{r) with a parameter at as 
a prefactor - summation is extended over the N elec- 
trons of the system. We assume that electrons with op- 
posite spin, but sharing the same spatial orbitals, have 
the same a^. In the cases where p-orbitals are needed we 
only use one parameter for all them, in order to preserve 
the spherical symmetry. For our light atoms we then have 

where nis(r), n2s{r) and n2p{r) arc the orbital densities, 
obtained squaring the appropriate atomic orbitals. All 
quantities only depend on the radial distance r. We have 
obtained the same qualitative behavior for the values of a 
as those obtained with the previous method and we can- 
not extract any additional information about the quality 
of the functionals in this way. As expected, we conclude 
that the use of the laplacian of the electron density is the 
correct way to generate a representative set of KEDs. 

For the sake of completion, we also present the re- 
sults obtained with a more sophisticated fully nonlocal 
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H 


He 




Li 


Be 




B 




c 




N 




o 




F 


Ne 


Ave. 


Range 


TF 





163 





163 





166 





169 





169 





169 





169 





169 





169 





168 


0.167 


0.006 


GEA2 





144 





144 





147 





149 





150 





151 





151 





152 





152 





152 


0.149 


0.008 


TF5W 





129 





129 





132 





135 





136 





138 





138 





140 





140 





140 


0.136 


0.011 


TFvW 





000 





000 





046 





066 





074 





080 





083 





089 





095 





097 


0.063 


0.097 


TF9W 





143 





143 





146 





148 





149 





150 





150 





151 





151 





151 


0.148 


0.008 


TF-N 





162 





162 





165 





167 





167 





168 





167 





167 





167 





167 


0.166 


0.006 


Pear 





164 





164 





168 





170 





171 





171 





171 





170 





170 





170 


0.169 


0.007 


DK 





144 





144 





149 





160 





159 





159 





158 





158 





157 





157 


0.155 


0.016 


LLP 





150 





150 





151 





153 





154 





155 





155 





155 





155 





155 


0.153 


0.006 


OLl 





143 





143 





146 





148 





149 





150 





151 





151 





151 





151 


0.148 


0.008 


OL2 





144 





144 





147 





149 





150 





151 





151 





152 





152 





152 


0.149 


0.008 


Thak 





148 





148 





150 





152 





153 





153 





153 





153 





154 





153 


0.152 


0.005 


B86A 





149 





149 





151 





153 





154 





154 





154 





154 





155 





154 


0.153 


0.005 


B86B 





147 





147 





148 





151 





151 





152 





152 





152 





152 





152 


0.150 


0.005 


DK87 





148 





148 





150 





152 





152 





153 





153 





152 





153 





152 


0.151 


0.005 


PW86 





150 





150 





152 





154 





155 





155 





155 





155 





155 





156 


0.154 


0.005 


FW91 





147 





147 





148 





151 





152 





152 





152 





152 





153 





152 


0.151 


0.006 


LG94 





148 





148 





149 





151 





152 





152 





152 





152 





152 





152 


0.151 


0.005 


vW 





000 





000 





000 





000 





002 





002 





002 





003 





003 





003 


0.002 


0.003 


ABSP 





000 





000 





004 





016 





025 





033 





040 





045 





049 





053 


0.026 


0.054 


GR 





000 





000 





006 





015 





023 





031 





037 





043 





048 





051 


0.025 


0.051 



TABLE III: Values of a when the semilocal functionals are fitted to tj^{r), Eq. (|17p . 







H 


He 




Li 


Be 




B 




C 




N 




O 




F 


Ne 


Ave. 


Range 


TF 





166 





166 





166 





168 





171 





167 





162 





161 





160 





160 





165 





012 


GEA2 





187 





187 





188 





190 





192 





188 





187 





186 





185 





186 





188 





007 


TF5W 





225 





225 





227 





228 





228 





225 





224 





221 





220 





221 





224 





009 


TFvW 





918 





918 





893 





873 





838 





805 





779 





750 





727 





710 





821 





208 


TF9W 





189 





189 





191 





193 





195 





190 





190 





188 





187 





188 





190 





007 


TF-N 





210 





213 





210 





207 





204 





197 





196 





191 





190 





191 





201 





023 


Pear 





166 





166 





166 





168 





170 





165 





162 





161 





160 





161 





164 





010 


DK 





218 





218 





207 





208 





209 





203 





200 





197 





194 





194 





205 





024 


LLP 





199 





199 





199 





199 





200 





195 





193 





191 





190 





191 





196 





010 


OLl 





192 





192 





193 





195 





197 





192 





191 





190 





189 





190 





192 





008 


OL2 





190 





190 





191 





194 





195 





191 





190 





188 





188 





189 





191 





008 


Thak 





203 





203 





202 





202 





202 





197 





195 





193 





191 





192 





198 





Oil 


B86A 





203 





203 





202 





202 





202 





196 





195 





192 





191 





192 





198 





Oil 


B86B 





217 





217 





216 





215 





215 





210 





208 





205 





204 





204 





211 





013 


DK87 





210 





210 





209 





208 





208 





203 





201 





198 





196 





197 





204 





013 


PW86 





201 





201 





200 





200 





200 





195 





193 





191 





190 





191 





196 





Oil 


PW91 





200 





200 





200 





200 





201 





196 





194 





192 





191 





191 





197 





010 


LG94 





207 





207 





207 





206 





206 





200 





198 





195 





193 





194 





201 





014 


vW 





000 





000 





027 





064 





104 





143 





179 





215 





247 





275 





125 





275 


ABSP 





378 





111 





034 





118 





173 





209 





232 





256 





271 





280 





206 





345 


GR 





252 





000 





044 





107 





156 





192 





217 





243 





260 





270 





174 





270 



TABLE IV: Values of a when the semilocal functionals are fitted to is'(r), Eq. (|17p . 



functional. By "fully nonlocal" we mean the functional 
explores the whole space when evaluating the contribu- 
tion to the kinetic energy from any point of the sys- 
tem. We have chosen the functional developed in 1985 
by Chacon- Alvarellos-Tarazon^, the simplest of a fam- 
ily of functionals.— '^"i^^'^^i^^ In Table fVl we show the rel- 
ative errors for the total kinetic energies and the values 
of a obtained with the aforementioned procedure and 
this functional approximation. Note the errors are about 
those obtained for most of the GGA functionals (in the 



average, smaller than 3.5%), and the values of a do not 
represent any clear improvement over the GGA ones. A 
complete study of the rest of the related fully nonlocal 
functionals, as well as a number of another kinetic en- 
ergy functionals, has been done in Ref. [6^ and will be 
presented elsewhere.— 
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H He Li Be B C N O F Ne 


Ave 


Range 


Rel. error 
a 


0.107 -0.041 -0.012 0.022 0.021 0.032 0.037 0.032 0.023 0.011 
0.107 0.055 0.082 0.105 0.122 0.138 0.155 0.170 0.186 0.205 


0.034 
0.133 


0.148 
0.150 



TABLE V: Relative errors in the total kinetic energy and values of o when adjusting with Eq. (|17|l for the original 1985 CAT 
functional. The average of the relative errors is made over their absolute values. 



VIII. A BRIEF GRAPHICAL STUDY 

Up to now, integrated values for the study of the KED 
have been discussed. Novif we present a graphical study 
of the KED for the neon atom; our aim is to show how 
the new technique developed works and the qualitative 
behavior of the approximate KEDs. 

For the sake of briefness, we have chosen three repre- 
sentative functionals: the TF functional, the GEA2 ap- 
proximation and the vW functional. TF yields the low- 
est values of the quality factor, GEA2 is representative of 
those GGA functionals that give errors of about 1% for 
the total kinetic energy and the lowest values of a of the 
usual GGA functionals. Finally, vW is also studied due 
to its theoretical importance and its special behavior. 

In Figs. [U [5] and [3] we show the orbital-based KED for 
the Ne atom as a thick solid line, the approximated KEDs 
corresponding to the three functionals as dashed lines and 
the approximate KED that includes the contribution due 
to the laplacian of the electron density is depicted with 
a thin solid line. 







4n\(r) 




' ' \ 
' ' \ 






h 













T ' 1 ' 1 1 1 1 1 ' 

0.0 0.2 0.4 0.6 0.8 1.0 

r (a.u.) 




-100- 

-200 -|U 1 1 1 1 1 1 1 1 1 

0.0 0.2 0.4 0.6 0.8 1.0 

r (a.u.) 

FIG. 1: Radial TF kinetic energy densities for the neon atom. 

Being t^^i^) always positive, we choose it for con- 
venience as the reference KED, whereas the correc- 



FIG. 2: Radial GEA2 kinetic energy densities for the neon 
atom. 



tions with the laplacians of the density are included 
as "laplacian contributions" to the approximated func- 
tionals. This can be done because comparing tj^{Y) = 
tg{r) + ato{r) with ^^""^(r) is equivalent to compare 
i^""'^(r) — afo(r) with tg{r). We use tg{r), always pos- 
itive, because it shows more clearly the KED in differ- 
ent regions of the space. Specifically, the figures exhibit 
a clear first shell and a small shoulder at the second 
shell (these features are verified against the radial den- 
sity, 47rr^n(r), and do not appear at the same positions 
when using tg^{r)). 

For the TF functional (Fig. [T]) we see that, without 
the laplacian contribution, the KED distribution has its 
main contribution to the kinetic energy in a region nearer 
to the nuclei than tg{r) has. The approximate KED with 
the laplacian contribution yields a more similar distribu- 
tion but two eye-catching pathologies are shown. The 
peak corresponding to the Is shell is exaggerated and a 
defect of KED is found near the nucleus. Instead, the 
asymptotic behavior seems to be almost correct. 

For GEA2 we can observe (Fig. [5]) qualitatively almost 
the same behavior than the TF, but now the pathologies 
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600 




-100- 

-200-1 , 1 , 1 , 1 , 1 , 

0.0 0.2 0.4 0.6 0.8 1.0 

r (a.u.) 



FIG. 3; Radial vW kinetic energy densities for the neon atom. 



are bigger. It is now clear how this gradient correction 
yields better values of the total kinetic energy than TF, 
without improving the local behavior of the KED, being 
the differences with the reference distribution bigger and 
the value of a worse. 

For the von Weizsacker functional (Fig. [3]) very dif- 
ferent results are obtained. The approximated KED is 
always smaller than the exact one, reflecting that when 
the von Weizsacker functional is written in the form given 
by the equation ([5]), it is a local lower bound for the first 
definition of the KED. That is, the vW functional is not 
only a lower bound to the total kinetic energy but is also 
a local lower bound for the KED;^ and there is no need 
for any laplacian correction. Note that in this case the 
contribution of the laplacian of the electron density yields 
a curve that is almost indistinguishable from the curve 
without it, as commented in section IVIIl and it is not 
depicted. It can be noted that tyw is only correct near 
the nucleus and in the asymptotic decay - the single or- 
bital regions-, as expected from the definition of the vW 
functional i2£ 



IX. CONCLUSIONS 

We have developed a method to test kinetic energy 
density functionals attending not only to the total kinetic 



energies but looking at the local behavior of the kinetic 
energy densities associated to the functionals. To obtain 
quantitative measure of the accumulated difference be- 
tween the distributions we have defined a quantity cr that 
we have called quality factor. Due to the non-uniqueness 
of the KED definition we have performed comparisons 
with an infinite family of kinetic energy densities, tg{r), 
generated by adding to the definition of the orbital-based 
KED the laplacian of the electron density multiplied by 
a variable prefactor. 

The procedure have been employed to test the local 
quality of twenty one semilocal functionals. We check 
the fitting of the KED corresponding to each functional 
to the closest tg{r) given by Eq. (fT7|) . For a given func- 
tional, we got values of a (that weights the contribution 
of the laplacian), corresponding to a minimum value of 
cr, that show a small dependence on the atomic num- 
ber Z. And for that value of a the corresponding KED 
is always closer to the mean of the first and second defi- 
nitions of the orbital-based KED than to the definitions 
themselves. This result recalls that this mean has been 

proved to be the more natural definition for a "classical" 
KED.20-63 

The main result found is the unexpected failure of all 
semilocal functionals but those with a full vW term to 
improve the local pathologies of the Thomas-Fermi func- 
tional. Our measurement technique assures that, even 
the GGA corrections to the TF functional always yield 
better total kinetic energies, this semilocal functionals 
get worse local KEDs. This result confirms the prelimi- 
nar calculations presented in Ref . [gsI. The only exception 
to this rule is the Pearson functional, that yields bad val- 
ues for the total kinetic energies. 

Finally, we have qualitatively studied the behavior of 
the KEDs, showing the characteristic pathologies of the 
TF functional, that exhibits an excess in the first peak 
of the density (corresponding to the Is orbital) and a 
defect near the nucleus. This pathologies are always en- 
larged in all the semilocal functionals, reflecting the main 
conclusion of this paper. 
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